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We theoretically investigate a generalized "which-path" measurement on an electronic Mach- 
Zehnder Interferometer (MZI) implemented via Coulomb coupling to a second electronic MZI act- 
ing as a detector. The use of contextual values, or generalized eigenvalues, enables the precise 
construction of which-path operator averages that are valid for any measurement strength from the 
available drain currents. The form of the contextual values provides direct physical insight about 
the measurement being performed, providing information about the correlation strength between 
system and detector, the measurement inefficiency, and the proper background removal. We find 
that the detector interferometer must display maximal wave-like behavior to optimally measure 
the particle-like which-path information in the system interferometer, demonstrating wave-particle 
complementarity between the system and detector. We also find that the degree of quantum erasure 
that can be achieved by conditioning on a specific detector drain is directly related to the ambiguity 
of the measurement. Finally, conditioning the which-path averages on a particular system drain 
using the zero frequency cross-correlations produces conditioned averages that can become anoma- 
lously large due to quantum interference; the weak coupling limit of these conditioned averages can 
produce both weak values and detector-dependent semi-weak values. 
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I. INTRODUCTION 

The construction of electronic Mach-Zehnder interfer- 
ometers (MZIs) in the solid state is a recent innovation 
in the fabrication and control of coherent mesoscopic sys- 
tems. The first experiment of this kind, published by 
the Heiblum group^ used the edge states^ of an integer 
quantum Hall Corbino geometry as the electronic analog 
of light beams and quantum point contacts (QPCs) 3 ~ 5 
as the electronic analogs of optical beam splitters to con- 
struct an interferometer with a visibility of 62%. Other 
interferometer designs have since been similarly imple- 
mented as electronic interferometers in the integer quan- 
tum Hall regime^— 

The electronic MZI differs from its optical counterpart 
in several respects. The arms of the MZI accumulate 
relative phase differences not only due to kinetic prop- 
agation of electrons along the arms, but also because 
the electrons are charged particles and can thus acquire 
a geometric Aharanov-Bohm phas o 10 ' 11 when the arms 
enclose a magnetic flux. This charge can also lead to 
strong electron-electron interactions, giving rise to a vari- 
ety of effects that have no counterpart in an optical MZI. 
For example, the interactions can produce differences 
in the counting statistics ) 12 ' 13 can induce temperature- 
dependent decoherencej^— can be used to detect exter- 
nal charges ; 20 ' 21 and can even lead to lobe structure in 
the visibility at high voltage biasj2£~— Such lobe struc- 
ture was unexpected and has generated numerous the- 
oretical explanations some hypothesizing Luttinger 
liquid physics as the cause. 

Here we take a more modest theoretical approach for 
describing electronic MZIs in a quantum Hall system that 
focuses on the low bias regime within a single-particle 
edge-state model. We consider a configuration of two 
such single-particle electronic MZIs coupled together by 



the Coulomb interaction, as shown in Figure [T] This ge- 
ometry has similarities to Hardy's parados 3 ^— and has 
been considered previously at various levels of detail by 
several authors 3 ^— to explore such phenomena as quan- 
tum erasure^ 2 - - — and Bell inequality violations. In our 
treatment, we include interactions between the MZIs via 
a minimal phenomenological model that adds a relative 
interaction-induced phase shift between a pair of elec- 
trons that occupy adjacent edge states simultaneously. 
The relative phase shift has the effect of entangling^ - — 
the states, mixing the path information of the two MZIs. 




FIG. 1. (color online) Schematic of coupled electronic MZIs. 
An ohmic source in a quantum Hall system at filling factor 
v = 2 injects chiral excitation pairs biased at energy eV rel- 
ative to the ohmic reference drains Di, D2, Si, and 5*2 into 
independent edge channels. The bias is kept low enough to 
allow only one excitation per channel on average. The outer 
(red) channel is transmitted entirely through the QPC Qf and 
partially transmitted through Qf and Qf, forming the system 
MZI. The inner (blue) channel is reflected entirely from Q\ 
and Q2 and partially transmitted through Qf and Q%, forming 
a separate detector MZI. The Coulomb interaction between 
the copropagating arms L d and U s induces an average rela- 
tive phase shift 7 between each excitation pair that couples 
the interferometers. 
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This kind of interaction has been experimentally shown 
to be capable of producing a 7r-phase shift on a single 
electron — perhaps the most dramatic difference from the 
optical analog, 1 - Thus, all elements of our theoretical 
analysis are based on currently available technology. 

Our work considers the task of detecting which-path 
information in one MZI by using the second MZI as the 
detector. Since the system and the detector are identi- 
cal devices, this arrangement has several appealing fea- 
tures. First, the symmetry of the geometry indicates 
that there should be a duality between "which-path" 
information in one MZI versus "which-fringe" informa- 
tion in the other. We will show that this is indeed 
the case, which relates this work to earlier "controlled 
dephasing" experiments We apply the contextual 
values formalis m 55 ' 56 for generalized measurements 5 -^— 
to show how even with inefficient detectio n 60 ' 61 and low 
visibility, the which-path information may be extracted 
from the detector currents systematically. Next, the fact 
that both the system and the detector have their own 
inputs, outputs and coherence allows the effects of mea- 
surement to be explored in detail. In particular, the cor- 
relations between them can be experimentally measured 
and analyzed, taking various forms such as joint counting 
statistics, or even conditioned measurements (see Sukho- 
rukov et al*£ for an example with incoherent electrons). 
The ability to condition (or post-select) measurements 
performed on a system with quantum coherence also 
allows the possibility of measuring weak valuesi 55 ' 63 " — 
Weak values, in addition to being of interest in their 
own right, have been shown to be useful as an amplifica- 
tion technique for measuring small variations of a system 
parameter— ~—, as well as for tests of bona fide quantum 
behaviorS^-Ii 

The paper is organized as follows. In Section[Hl we de- 
scribe the considered geometry and our physical model- 
ing principles before reviewing and applying the scatter- 
ing theory of mesoscopic transport to the two edge state 
physics. In Section IIII1 we take the joint predictions of 
the scattering model and interpret them as a measure- 
ment by the detector electron that extracts information 
about the path state of the system electron. We describe 
a principled way to construct this information from the 
data in a variety of regimes. In Section IIVI we intro- 
duce the conditioning procedure and use it to clarify the 
phenomenon of quantum erasure, as well as to calculate 
conditioned averages of the which-path measurements. 
We shall see that both weak values and semi- weak values 
will compete as weak coupling limits of these conditioned 
averages. In Section [V] we describe our conclusions. 



II. COUPLED MZIS 

We consider a pair of electronic MZIs embedded in 
a two-dimensional electron gas in the integer quantum 
Hall regime at filling factor v = 2 as illustrated in Fig- 
ure [1] An ohmic source with a small DC-bias eV above 



the Fermi energy Ep injects chiral electron-like edge ex- 
citations of charge e into the sample that propagate uni- 
directionally along two independent edge channels. 2 Each 
edge channel forms an interferometer from two appropri- 
ately tuned quantum point contacts (QPCs) that coher- 
ently split and then recombine the possible paths. The 
relative phase between the arms of each interferometer 
is determined not only by a local dynamical phase accu- 
mulated during kinetic propagation along each arm, but 
also by a global geometric phased (in the form of the 
Aharonov-Bohm (AB) effect arising from the closed 
paths. After the paths interfere, the charges are col- 
lected at ohmic reference drains held at the Fermi en- 
ergy, producing fluctuating output currents that can be 
temporally averaged. 

The two interferometers accrue an additional relative 
phase shift due to the Coulomb interaction where the 
charges copropagate. Intuitively, the mutual repulsion 
affects the dynamical phases by effectively warping the 
propagation paths, which also affects the geometrical 
phases by changing the areas enclosed by the paths. A 
more careful analysis of the joint interaction phase is pro- 
vided in the Appendix. Such additional relative phase 
has the effect of entangling the joint state of the two inter- 
ferometers, mixing the which-path information. Due to 
the entanglement, extracting information from the drains 
of one interferometer allows one to infer correlated which- 
path information about the other interferometer. That is, 
one interferometer can be used as a detector to indirectly 
measu re 55 ' 57 ~ — the which-path information of the other. 
As we shall see, the characteristics of the measurement 
will depend on the tuning of the detector interferometer 
as well as on the coupling phase. 

We model the coupled MZI system using the elastic 
scattering approach of Landauer and Biittiker— ~— for co- 
herent charge transport. As the transport is largely bal- 
listic in the integer quantum Hall regime, the formalism 
directly relates the average currents collected at each 
ohmic lead I £ {D±, D2, Si, S2} to the transmission prob- 
abilities Pi(E) £ [0, 1] for plane waves of fixed energy E 
to traverse the sample successfully. Treating the ohmic 
leads as thermal reservoirs, the average currents from 
spinless single-channel transport are, 

pe-V irp 

h=ej —(f(E + eV)-f(E))P l (E) (la) 
e 2 V 

*—Pi(E F ). (lb) 

Here, f(E) = (exp {{E - E F )/k B T) + l)" 1 is the equi- 
librium Fermi distribution relative to the Fermi energy 
Ep at a temperature T; h is Planck's constant; and, ks 
is Boltzmann's constant. 

The approximate equality (jlbl) holds in the low-bias 
regime when Ep 3> eV 3> ksT and the transmission 
probabilities Pi are constant with respect to the small 
variations in energy. We assume that the source oper- 
ates in such a regime. Due to the small spectral width 
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of the source, the fermionic excitations will then be well 
approximated as plane waves at a fixed energy on the 
scale of the sample; hence, on average only one excita- 
tion per channel will occupy the sample and intrachannel 
interactions can be ignored. In particular, we avoid the 
anomalous lobe structure in the interference that appears 
at higher bias i^2r— 

We also assume for simplicity of discussion that the 
source only injects spinless excitation pairs with one ex- 
citation per channel so that the coupling interaction be- 
tween the channels is fixed; the results will be averaged 
over a more realistic source distribution in Section UlI Gl 
With these approximations, the initial joint scattering 
state for an excitation pair can be written in second- 
quantized notation as, 




FIG. 2. (color online) Complementary QPC balance parame- 
ters ((4} for m £ {d, s} and i £ {1, 2} as parametrized by the 
balance angle 6™ . Left: transmission probability T™ (solid, 
blue) and reflection probability R™ (dashed, red). Right: 
particle-like parameter <5f (solid, blue) and wave-like param- 
eter e*™ (dashed, red). 



|*) = a dt a st |0), (2) 

where |0) is the filled Fermi sea of the edge channels and 
a rf t anc i g^t are creation operators for plane waves of a 
fixed energy injected into the inner and outer channels, 
respectively. Operators corresponding to different edge 
channels commute due to the independence of the chan- 
nels. 

The inner channel will form a Mach-Zehnder interfer- 
ometer as shown in Figure [TJ which we refer to as the up- 
per MZI, or the detector MZI. Similarly, the outer chan- 
nel will form an identical interferometer, which we refer 
to as the lower MZI, or the system MZI. We will use the 
lowercase superscripts d and s throughout to distinguish 
quantities specific to the detector and the system, re- 
spectively, and to avoid confusion with the detector and 
system drains that we denote with capital letters D\ , D2 , 
S\ and 52. 

The QPCs Qf, Qi, Qf, and Q| shown in Figure [Q 
each elastically scatter the plane waves, affecting only the 
complex amplitudes of the joint scattering state. Hence, 
for m £ {d, s}, i £ {1,2} we can represent the effect of 
each QPC as a unitary scattering matrix, 

j-j-rn _ I c L i ' i \ lo\ 

\e ix ^r™ e^?t™ ) ' K ' 

where t™ = y/T™ and r" 1 = i^R" 1 are given in terms of 
the transmission and reflection probabilities T™ £ [0, 1] 
and Rf = 1-T™ though Qf . The additional scattering 
phases xT and £f may arise from QPC asymmetry. 

The QPCs are kept tunable subject to the constraints 
that the outer channel is fully transmitted through Qf 
and Q2 and the inner channel is fully reflected from 
Qf and Q| to create the two separate interfering paths. 
There is an additional QPC near drain S± not shown in 
Figure Q] that is kept fixed to allow full transmission of 
the outer channel and full reflection of the inner channel 
in order to divert the outer channel for collection at the 
drain S±. 

For later convenience we also introduce the comple- 



mentary QPC balance parameters, 

<*r = 7T-iCe[-M], (4a) 
ef = 2^Tf r W£ [0,1], (4b) 

for m £ {d, s} and i £ {1, 2} that satisfy (ef ) 2 + (<5f ) 2 = 
1. All such QPC parameters can be related by a QPC 
balance angle 0f £ [0,tt/2] such that, If 1 = cos 2 6»f, 
Rf 1 = sin 2 !?!", 5™ = cos 20™, and e™ = | sin 2(9™ | as il- 
lustrated in Figure [2] We shall see that the parameters 
5™ control the particle-like path-bias of the excitation 
after a QPC, while the parameters e™ control the com- 
plementary wave-like interference visibility. 

The joint state ([2]) can be scattered through Qf and 
Ql using into the basis of the MZI paths, yielding the 
replacements, 

a^=Mt d ia\ d +e^rUU ( 5a ) 
a st = e lx h s ia [ s +e z ^r s 1 al s . (5b) 

During propagation to the second pair of QPCs, each 
path p £ {L d , U d , L s , II s } accumulates an additional dy- 
namical phase <j) p that depends on the excitation energy 
and the path-length. When the paths recombine, the dif- 
ference between the dynamical phases contributes to the 
interference. Closing the paths for MZI m £ {s, d} also 
contributes a relative geometric Aharonov-Bohm (AB) 
phase 4>™b that depends on the magnetic flux enclosed 
by the path. 

We can compactly account for the various phase effects 
contributing to the interference by defining tuning phases 
for each MZI, 

d = 4> d AB + <j> Ld - 4> UA + xi - £1 , (6a) 
s - f AB + LS - fo. +xl- Ci ■ (6b) 

Finally, the joint scattering amplitude correspond- 
ing to co-occupation of L d and U 8 acquires a effective 
Coulomb interaction phase 7 that couples the two in- 
terferometers. (See the Appendix for discussion about 
how the Coulomb effect can produce such a phase shift.) 
This interaction phase compactly encodes opposing shifts 
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in the combined dynamical and geometric phases of each 
MZI due to the Coulomb repulsion of the charge pair. 
For simplicity, we assume for now that the relative phase 
is constant; we will allow it to fluctuate for consecu- 
tive pairs in Section IIII Gl We also note that any ad- 
ditional Coulomb phase acquired during copropagation 
after QPC Qf and before QPC Q\ will only contribute 
to the tuning phase <fi d and can therefore be ignored. 

After adding the phenomenological phases, the scat- 
tered joint state just before the second pair of QPCs is, 

|*') = (tft{ e^ d +^a{ d al +r( rj^o^ (7) 

+ rf t{ e^a f ud al + tf rf e^+^a^a^) |0), 

up to a global phase of exp(i(<fo/d -I- (j>u' + S,f + £!)) not 
written. 

The interaction phase 7 has the effect of entangling the 
two interferometers, which we can show by computing the 
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(8) 



We see that the entanglement reaches a maximum when 
the phase 7 — > it and vanishes as 7 — > 0. Further- 
more, the entanglement directly depends on the QPCs 
Qf and Qf preparing interfering wave-like excitations 
in each MZI, which is measured by ef ef ; maximum 
entanglement can only occur for balanced QPCs with 
Tf = T{ = 1/2, or efef = 1. 

At this point, we conceptually break the symmetry be- 
tween the two interferometers to treat one as a detector 
for information about the other. We will treat the upper 
MZI as the detector and the lower MZI as the system be- 
ing measured, though obviously we could exchange those 
roles by the symmetry of the geometry. To do this we 
finish scattering the detector MZI through Q 2 into the 
basis of the ohmic detector drains {Di, D 2 } using ((3]), 



(9) 



l D 2 , 



(10) 



yielding, 

l*"> - («k 

+ cJ D2 [C D2 , L se l * s t{ a[ s + C D2 j Js r{ fit,.]) |0), 

up to the same global phase as in (JTJ). For later con- 
venience, we have defined the detector scattering ampli- 
tudes, 

C Dl , L . = M (tftie^" +r d x r() , (11a) 
C DuU . = e lxi (tftieW+rt + rfr d ^ , (lib) 



C D2 , L s=e^ (ttrie^+rfti), (11c) 
' (tfrU^^+rfti). (lid) 



C Da ,u> = e« 



III. MEASUREMENT INTERPRETATION 

The joint scattering model is useful for computing 
probabilities and average currents, but it does not pro- 
vide direct insight into the measurement being performed 
by one interferometer on the other. To make the con- 
nection to measurement more apparent, we will use the 
contextual values formalism 5 -' 56 that links the detector 
drain probabilities directly to the "which-path" operator 
for the system. We will see that we can understand the 
various subtleties of the measurement quite transparently 
using this technique. 



A. POVM 

To facilitate the interpretation of the distinguishable 
detector drains as the outcomes of a measurement being 
performed on the system, we define the single particle 
state kets, 

\Di) = 4)J0), |£ 2 )=4 2 |0), (12a) 

\L s )=al\0), \U s )=aU0), (12b) 

define the reduced system state in absence of interaction, 

\r) = e^t\\L s )+r{\U s ), (13) 

and write (|10l) in the form, 

|*"> = \Dt) ® M Dl \r) + \D 2 ) ® M D2 \tP s )- (14) 

The interaction with the detector in (fT4"|) is entirely 
represented by operators acting on the reduced system 
state (|13p that contain all the scattering information of 
the detector, 

M Dl =C Dl , L s\L s )(L s \ + C DuU s\U s )(U s \, (15a) 
M D2 =C D2 , L s\L s )(L s \ + C D2 ,us\U s )(U s \. (15b) 

The operator Mo 1 encodes the interaction followed by 
the absorption of the detector excitation at the drain 
D\. Similarly, the operator Md 2 encodes the interaction 
followed by absorption at D 2 - We refer to {Mb^MbJ 
as measurement operators*^— 

As the coupling phase 7—^0 the measurement op- 
erators (|15p become nearly proportional to the identity. 
We call 7^0 the weak coupling limit since the reduced 
system state is only weakly perturbed for small 7. Con- 
versely the limit 7 — > tt is called the strong coupling limit 
since the measurement operators are maximally different 
from the identity and maximally perturb the reduced sys- 
tem state. 

The measurement operators also form a Positive 
Operator- Valued Measure (POVM) on the system, 

E Dl = Mp i M Dl , E D2 = M ^ M D2 , (16) 
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such that E Dl + E D2 = 1. The POVM elements 
{Ed 11 Ed 2 } act as probability operators for the measure- 
ment outcomes. 

Hence, the probability of absorbing the detector exci- 
tation at a drain D g {£>i,£>2} can be expressed either 
as an expectation of the projection operator of the de- 
tector drain under the joint state ()10j) or, equivalently, 
as an expectation of the probability operator (|16p under 
the unperturbed system state (fT3l) . 



(17) 



p D = \(D\y"}\ 2 = (r\E D m 

= \C D ^t d \ 2 + \C D ^r d \ 2 . 



By working with the reduced state (TT31) . the measure- 
ment operators (fl"5j) . and the probability operators (fl6|) . 
we treat the detector as an abstract entity whose sole 
purpose is to measure the system. Such abstraction al- 
lows us to more clearly examine the measurement being 
made upon the system. 



B. Contextual Values 

In order to relate the measurement on the system to 
observable information that we can interpret, we will use 
contextual value s 55 ^ 56 to formally construct system ob- 
servables from the probability operators (|16p . This for- 
malism acknowledges that the only quantities to which 
we have experimental access are the detector drain prob- 
abilities, so all observations we wish to make about the 
system must be contained somehow in those probabili- 
ties. Generally, the correspondence between the detector 
drains and a particular system observable will be imper- 
fect, but we can compensate for such ambiguity of the 
detection by weighting the drain probabilities with ap- 
propriate values for the particular measurement setup. 

Generally, we cannot construct information about just 
any system observable from a particular measurement. 
To find which observables we can measure, it is useful 
to decompose the probability operators ([TB| into an or- 
thonormal basis for the observable space. In our case, the 
system state space is two-dimensional, so any Hcrmitian 
operator can be spanned by the four basis operators, 

a s = l= \L S )(L S \ + \U S )(U% (18a) 

d> 1 =dZ = \L'){U'\ + \U'){L'\, (18b) 

a s 2 =a s y = -i(\L s }(U s \-\U s )(L s \), (18c) 

a s 3 =a s z = \L*)(L s \-\U s }(U s \, (18d) 

which are equivalent to the identity operator and the 
Pauli spin operators. To find the real components of 
an observable in this basis we introduce the normalized 
Hilbert-Schmidt inner product between operators, 



A, B 



Tr A^B 



under which the operator basis is orthonormal, 

(3£,*;) = ^. (20) 

Here //, v € 0, 1, 2, 3 and 6^ u is the Kronecker delta that is 
1 if fx = v and otherwise. Using this basis, any system 
observable can be written, 



a lt = (A,df t )=Tt 



Aaf 



A 



(21a) 
(21b) 



where {a M } are real- valued components of the observable. 

Using (|2"Tj) . we can expand the probability operators 
on the system (TTH1) in the basis (fTE)) to determine their 
structure, 



E Dl =\{P%- V d A d ) a s - l -V d T d <7|, (22a) 
E D . 2 = \{pt + V d A d ) a s + ]-V d T d a s 3 , (22b) 



where we see that the measurement is characterized by 
the detector parameters, 

P% = 2(T?T$ +RtRi) =1 + 6(8$, (23a) 
Pt =2(T?R d +R d T d ) = l-5 d 5 d , (23b) 

V d =A^T d R d T d R d = 44, (23c) 

T d = sin ^ sin Q + (23d) 

A d = cos<?(> d -r d . (23e) 

defined in terms of the QPC balance parameters (HJ, the 
tuning phases ([6|), and the coupling phase 7. We will 
describe these parameters in detail in the next section. 

The probability operators only contain components in 
the subspace spanned by {ctq, ct|}; therefore, we can only 
construct observables that are contained within that sub- 
space. That is, we can construct any observable of the 
form A — aQ&Q+a^a^. We denote observables of this form 
as being compatible with the measurement (|22p . Other 
observables are incompatible with the measurement. 

To construct such a compatible system observable from 
the measurement, we expand its operator directly in 
terms of the probability operators (pi 



Tr[i] ' 



(19) 



A = a (To + 0,30-3 = a Dl E Dl + a D2 E D2 . (24) 

The required expansion coefficients ao 1 and ao 2 are gen- 
eralized eigenvalues, or contextual value s 55 ! 56 , of the op- 
erator. Using this expansion, we can recover the same in- 
formation on average as a projective measurement by us- 
ing only the drain probabilities, (A) = an 1 PD 1 +C(D 2 PD 2 - 
To determine the appropriate contextual values to as- 
sign in order to construct A, we insert (|22[) into (IM1) and 
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solve it as a standard matrix equation using the orthonor- 
mal basis, which yields the unique contextual values, 

* Dl = a - g + , (25a) 

a D2 = a + g - A d ) ■ (25b) 

As long as the contextual values do not diverge, the ex- 
pansion (IM|) of the compatible operator A is well defined, 
and we can perfectly recover its average, 

(A) = (r\A\r)=a + a 3 6{. (26) 

The observable parameter ao sets the reference point for 
the average, so contributes no information about the sys- 
tem; we will set it to zero in what follows without loss 
of generality. Similarly, the remaining parameter 03 sets 
the scale of the average; we will set it to one in what 
follows. 

We formally conclude that the detector drains perform 
a generalized measurement of the which-path operator <rj , 
as might be intuitively expected from the path-dependent 
interaction. Moreover, the QPC parameter Sf defined in 
(|4"aj) determines the particle-like which-path behavior on 
average. No other information about the system can be 
inferred from the measurement. 

The contextual values (|25|) are shown in Figure [3] for 
a few parameter choices. If they are equal to the eigen- 
values of the which-path operator, Od^Od-j = ±1, then 
the measurement is unambiguous: one obtains perfect 
knowledge about the path information with every drain 
detection, and the system state is projected to a pure 
path state. If the contextual values diverge, an 2 — >• 
±00, then the measurement is completely ambiguous: no 
knowledge about the path information can be obtained, 
and the system state is unprojected; however, will shall 
see in Section IIII Dl that the system state may still be 
unitarily perturbed by the coupling. In between these 
extremes the measurement is partially ambiguous: par- 
tial knowledge is obtained about the path information 
with each drain detection, and the system state is par- 
tially projected toward a particular path state. 

C. Parameters 

To better understand the parameters (|2"5|) we write the 
drain probabilities (fTTj) explicitly, 

P Dl = \ {Pi - V d ( A d + Sf T d ) ) , (27a) 

Pd 2 = \{PL + V d (A d + S[ T d )) . (27b) 

The probability Pd 1 is illustrated in Figured] for several 
values of the coupling strength. 

The particle-like parameters /8+,/?l € [0, 1] are deter- 
mined entirely by the path-bias parameters 8f and S d ; 
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FIG. 3. (color online) The contextual values (I25|l of the 
which-path operator 0-3: ctD 1 (solid, blue) and qd 2 (dashed, 
red), as a function of the coupling phase 7. The curves are 
shown for efficient detection V d = 1 and detector tunings 
4> d = {0, 7r/2, 37r/4, 7r}. The tuning strongly affects the ambi- 
guity of the measurement; moreover, the roles of the detector 
drains flip as the tuning varies from <f> d — to <f> d = 7r. 




FIG. 4. (color online) The drain probability Pd 1 P7)l as a 
function of the detector tuning <j> d and the which-path infor- 
mation <5f, shown for efficient detection V d = 1 and coupling 
phases 7 = {0, 7r/4, tt/2, it}. For zero coupling the interfer- 
ence is independent of the which-path information; for strong 
coupling 7 = 7r the interference maximally corresponds to the 
which-path information. 



they indicate the average background signal of each de- 
tector drain and satisfy (J3+ + j3 d _)/2 = 1. The wave- 
like parameter V d € [0,1] is determined entirely by the 
path-uncertainty parameters ef and ; it indicates the 
visibility of the interference. The parameter T d £ [— 1,1] 
indicates the deviation in the interference caused by the 



7 



coupling phase 7, which is the only effect of the charge 
coupling. The parameter A d 6 [—1,1] indicates the inter- 
ference unrelated to the path information of the system. 
As the coupling 7 — > 0, then T d — > and A d — » cosc/f> d , 
which recovers the signal for an isolated interferometer X 
As the coupling 7 — >■ it, then T d — > cos0 d and A d — > 0, 
and the interference maximally corresponds to the path 
information. 

The parameters (|2"5|) also give insight into the nature of 
the measurement by the role they play in the contextual 
values (|25|) . The parameter T d indicates the correlation 
between the detector drains and the which-path informa- 
tion. Its magnitude \T d \ e [0, 1] denotes the correlation 
strength, with 1 indicating perfect correlation and in- 
dicating no correlation; due to the inverse dependence in 
(|25p. any imperfect correlation will amplify the contex- 
tual values to compensate for the resulting measurement 
ambiguity. The sign of T d indicates the correspondence of 
the detector drains to the which path information, with 
— denoting the mapping {.Di,!^} ^ {L S ,U S } and + 
denoting the mapping {D2,Di} ■<-» {L S ,U S }. Note that 
the correlation strength depends not only on the coupling 
phase 7, but also on the tuning phase tfi d ; hence, it is pos- 
sible for the detector drains to be uncorrelated with the 
system paths even under strong coupling (e.g. examine 
cf> d = 7r/2 in Figure 2] when 7 = tt). 

The parameters f}+ and f3_ in (|25|) counterbalance the 
bias in the average drain background caused by a pre- 
ferred particle-like path. For instance, if ffi > f3_ then 
the signal at drain D\ is stronger on average in (|27[); 
hence, the contextual value (|25a[) assigned to D\ is pro- 
portional to the smaller value (3_ to compensate. 

The visibility parameter V d controls the wave-like in- 
terference produced by Qf and Q 4 -. The transmission of 
each QPC should be balanced in order to provide the in- 
teraction phase with an equal-amplitude reference phase 
for later interference. Any imbalance leads to inefficiency 
of the measurement^ by reducing the interference visibil- 
ity, effectively hiding the correlations. Such inefficiency 
increases the measurement ambiguity and results in an 
amplification of the contextual values. All correlations 
are hidden at zero interference visibility when either Qf 
or Q d is fully transmissive or reflective, T d ,T d 6 {0, 1}, 
which leads to divergent contextual values. Maximum in- 
terference visibility occurs for balanced transmission with 
V d = 1. We see that to optimally measure the particle- 
like which-path information for the system, the detector 
must itself exhibit maximal wave-like interference; the 
detector and system behaviors are therefore complemen- 
tary. 

The parameter A d is the portion of the interference not 
affected by the coupling, meaning A d + T d = cos<fi d . It 
indicates an additional bias in the drain correspondence 
caused by the interference not pertinent to the which- 
path measurement. The contextual values naturally sub- 
tract the contribution from this irrelevant background 
interference to retrieve the measurement information. In 
the limit of strong coupling 7 — > n, all the detector in- 



terference encodes the measurement result: T — > cos <j> 
and A d -> 0. 

Practically speaking, the detector must be calibrated 
in the laboratory before it can be used to probe an un- 
known system state. That is, the detector parameters 
(1231) must be predetermined by examining the drain out- 
puts of the detector under known system configurations. 
For example, pinching off QPC Qf to prevent any interac- 
tions allows most of the parameters to be set directly by 
tuning the detector QPCs and the magnetic field. The re- 
maining interaction parameter 7 can be inferred from an 
additional reference system state. Therefore, the process 
of detector calibration can be viewed as the experimen- 
tal determination of the appropriate contextual values to 
assign to the detection apparatus. 

D. Measurement Disturbance 

The measurement necessarily disturbs the system state 
by extracting information. We can see the effect of such 
disturbance by characterizing the system interferometer 
with analogous parameters to (f2"3"]). 





= 2 (T-f T£ +R S 1 R 2 , ) = 1 + 
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= cos ^ ~v s . 




(28e) 



Using these parameters the absorption probabilities for 
the system drain take the simple form similar to ([2Tj> . 

P Sl = \{P S + -V S {A S -8 d T s )), (29a) 
Ps 2 = \{P s - + V s (A s -5 d T s )). (29b) 

With efficient detection V d — 1 and strong coupling 
7 — > 7T, then 8f — > and A s — > 0, so the system drain 
probabilities display no interference, Pg 1 — > /3+/2 and 
Ps 2 — > /3£./2; that is, a strongly coupled, efficient which- 
path measurement will force particle-like statistics in the 
system pi&2£ 

The measurement disturbance may be analyzed more 
explicitly by rewriting the measurement operators (|15[) 
for the case of efficient detection V d = 1, 

M Dl = i e ix * e 1 ^' 2 \] 1 E)£ , (30a) 
M D2 = i e «i e ^ d / 2 U 1 E^l , (30b) 

^ = sm^ \ L '){L'\ + sin ^±2 \U')(U'\, (30c) 

E l £ = cos ^ I L s ) (L s I + cos ^-±1 \U*)(U S \, (30d) 
C/ 7 =exp(^|C/ s )(C/ s |) . (30e) 
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The disturbance manifests itself as two distinct processes. 
First, the positive roots of the POVM (p]) {E 1 ^, E 1 ^} 
perform the information extraction necessary for the 
measurement, partially projecting the reduced system 
state toward a particular path. Second, the coupling- 
dependent unitary factor U~ t contributes an additional 
evolution of the system that is unrelated to the extraction 
of information. The remaining phase factors contribute 
only to the global phase of the measured state and do 
not alter the subsequent measurement statistics. 

Unambiguous measurements extract maximal informa- 
tion from the system and thus project the system state 
to a definite path; they are frequently known as projec- 
tive or strong measurements. Ambiguous measurements 
extract partial information from the system and thus par- 
tially project the system state toward a particular path. 
Completely ambiguous measurements extract no infor- 
mation from the system and thus are completely unitary. 
When the system state is nearly unperturbed up to a 
global phase, the measurement is called weak, which cor- 
responds to the case of a nearly completely ambiguous 
measurement with a negligible unitary evolution. 



E. Strong Coupling 

An unambiguous measurement can only be obtained 
in the limits of efficient detection V d — > 1 and strong 
coupling 7 — > 7r. In this situation, the ambiguity will be 
determined only by the tuning phase of the detector (j) d , 
and the POVM will have the most symmetric dependence 
on the which-path operator, 



O-Di 

E Dl 

En.-, 



-1 



| (1-^3 W) 

l(l + al cos0 d ) 



(31a) 
(31b) 
(31c) 
(31d) 



As the tuning phase 4> d varies, the POVM elements 
oscillate between pure path projections and the identity, 
despite the strong coupling. The tuning-dependent drain 
ambiguity contributes to the inefficiency of the measure- 
ment by erasing the potentially extractable which-path 
information from the detector state. Indeed, we shall see 
in Section [IV Bl that such ambiguity in the measurement 
allows the system interference to be recovered by condi- 
tioning the system results on specific detector outcomes: 
Such a phenomenon is known as quantum erasure£2-~— 

In a laboratory quantum Hall system the AB phase will 
precess due to slow decay of the transverse magnetic field, 
so the tuning phase (f> d will also precess slowly. Hence, 
the ambiguity of the measurement will generally oscillate 
between extremes, while also flipping the correspondence 
of the drains to the which-path information. Despite any 



ambiguity in the measurement, however, the system will 
always be perturbed by the additional unitary evolution 
(|30ej) that induces a relative phase shift of ir/2 between 
the arms. Since the system state will be appreciably al- 
tered by the strong coupling, the measurement will not 
be weak even when completely ambiguous. Hence, ambi- 
guity of the measurement need not indicate weakness of 
the measurement. 

The measurement becomes unambiguous when the 
tuning is held fixed at cos(j) d = ±1. In this situation, the 
detector drains are perfect "bright" and "dark" ports: 
detection at the dark port will occur deterministically 
when the system excitation is in the upper arm. The de- 
tector drains are perfectly correlated to the which-path 
information of the system, so the system state is pro- 
jected to a definite path, and the measurement is strong. 



F. Weak Coupling Limit 

The weak coupling limit is the limit as the coupling 
phase 7 — > and the system and detector become nearly 
uncoupled. Since at zero coupling the measurement op- 
erators pOj) must either be zero or be proportional to the 
identity, the weak coupling limit of a measurement must 
have outcomes that are inherently ambiguous. Hence, we 
expect the contextual values (|25|) to diverge. However, 
since T d = sin (j) d {^/2) + cos0 d (7/2) 2 + 0(7 3 ) the nature 
of the divergence will also depend upon the tuning. 

If the tuning is not an integer multiple of tt, then both 
measurement operators fl30p will approach the identity 
as 7 — > and the measurement will be weak for all out- 
comes. That is, the system state will be nearly unper- 
turbed for any outcome of the measurement. In this case, 
T d = sin (f> d (7 / '2) + 0(7 2 ) and the divergence of the con- 
textual values (|25|) will be linear in 7. For an efficient 
detector with V d — 1, we find to 0(7 2 ), 



a Dl -> 1 



2 1 + cos d> d 



7 sin 0" 

2 1 - cos 4> d 
an 2 — > 1 H 



7 sin < 



Er 



I^i + 2sin^)<tn 

1 + C0S ^i-?Wi 



E D2 - ' y- l-^smflU'W 



(32a) 
(32b) 
(32c) 
(32d) 



The POVM has simple dependence on the projection to 
the upper path, which can also be written in terms of the 
which-path operator, \U S )(U S \ = (1 — 0"|)/2. The most 
symmetric case of (f> d = ir/2 is shown in the upper-right 
of Figure [3j 

However, if the tuning <fi d = nir with integer n, then 
only one of the measurement operators will approach the 
identity as 7 — > 0. The remaining outcome remains pro- 
portional to a projector with a vanishing coefficient and 
will thus strongly perturb the system state. Hence, only 
one contextual value diverges while the other remains a 
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constant eigenvalue. In this case, T d = (— 1)™ sin 2 (7/2) 
so the divergence will be quadratic in 7. We call such 
a measurement a semi-weak measuremenl^ since only a 
subset of outcomes are weak. For an efficient detector 
with V d = 1, we find, 

a Dl -> —i- ((-I)" + cos 2 2) , (33a) 
sm ^ \ 2, J 

a D2 -> -1- ((-!)» - cos 2 J) , (33b) 
sm ^ V 2/ 

E Dl -> ^(1 - (-1)™)! + (-I)"sin 2 ||£/')(^|, (33c) 
^ 5(1 + (-!)")! - (-l)"sin 2 1\U'){U'\. (33d) 

The POVM retains the simple dependence on the projec- 
tion to the upper path. The cases for n — and n = 1 
are shown in the upper-left and lower-right of Figure EH 
respectively. 

For the semi- weak measurement, the effect of absorp- 
tion at one of the drains is projective. The projective 
drain outcome unambiguously indicates that the system 
excitation took the upper II s path; therefore, the con- 
textual value assigned to the complementary drain is 
an eigenvalue. In contrast, the effect of absorption at 
the complementary drain only weakly perturbs the sys- 
tem state. Its outcome only ambiguously corresponds to 
which-path information; therefore, the contextual value 
assigned to the projective drain must be amplified. Such 
complementary behavior of the contextual value ampli- 
fication can be counter-intuitive, but it emphasizes that 
the function of the amplification is to compensate for the 
ambiguity of the measurement. 

We shall see in Section H V F I that while conditioned av- 
erages of the weak measurements ([32)) will lead to weak 
values, the conditioned averages of the semi-weak mea- 
surements (|3"3"|) have different limiting behavior and lead 
to different values. The two limiting cases will compete 
depending on the relative magnitudes of 7 and <fi d . 

G. Fluctuating Coupling 

If the coupling between excitation pairs is not a con- 
stant relative phase 7, but instead can fluctuate within a 
finite uncertainty width a around an average 7, then the 
average measurement will be correspondingly more am- 
biguous. We could quantify this effect by averaging the 
joint state (fit)]) over a range of coupling phases to create a 
mixed state represented by a density operator; the mea- 
surement operators (|15p and resulting POVM (TT6T) could 
then be generalized to an averaged measurement from 
that density operator. However, that procedure would 
be completely equivalent to the simpler procedure of av- 
eraging the probability operators (fTB"|) over the coupling 
width directly, which we choose to do here. 

For simplicity, we consider as a coupling distribution 
the raised cosine distribution, which is Gaussian-like, but 



P(y') l(fr) 




FIG. 5. (color online) Left: The raised cosine distribution 
showing a spread in the coupling phase centered at 7 = n/2 
by a half-width a — n/4. Right: The inefficiency factor n(a) 
defined in (|36[) as a function of the half-width a. 

has compact support. We center the distribution around 
7 G [0, 2n], and give it the half- width a G [0,7r]. The 
density for the distribution is nonzero in the domain 7' € 
[7 — a, 7 + a] and has the form, 

W) = ^(l + cos (-(-/- 7))). (34) 

An example of the distribution is shown in Figure [5] cen- 
tered at 7 — it/2 and with half- width a — tt/4. 

Averaging the probability operators (| 16[) only affects 
the constant T d , which is replaced by the averaged ver- 
sion, 

T d (7) f"^ di P{i)Y\i) = 77(^(7), (35) 

J '■f — cr 

( k 2 \ ( sincr\ , . 

•rf'H^H— )■ (36) 

A plot of the damping factor 77(17) is shown in Figure [SJ 
We assumed in (|2|) that the source emits only excita- 
tion pairs. However, any contribution of unpaired exci- 
tations in the initial joint state is equivalent to a contri- 
bution of joint states with 7 = 0. The net effect of the 
source emitting such unpaired excitations is thus to mod- 
ify T d by an additional probability factor P p that denotes 
the likelihood of pair emission. Hence, the only effect of 
an imperfect source is to introduce a net inefficiency fac- 
tor 77' = P p rj(a) € [0, 1] in the parameter T d . 

Since the contextual values (|25|) inversely depend on 
T d , any such inefficiency will introduce an overall ampli- 
fication factor of 1/77'. In other words, any uncertainty 
in the coupling strength will lead to additional ambiguity 
in the average measurement by degrading the portion of 
the detector interference that is coupled to the system. 

H. Observation Time 

Since ambiguous measurements provide less informa- 
tion about an observable per measurement, more mea- 
surements will be required to achieve a desired precision 
for an observable average. We can characterize the nec- 
essary increase in observation time as follows. The total 
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detector current is / = (e 2 VP)/h = e/r m according to 
(|lb[) . where P is the total probability for excitations to 
traverse the sample; hence, we can infer that the aver- 
age time per detector absorption is r m = h/eVP. For 
our single- particle model to apply we wish for the volt- 
age bias V to be low enough that the interferometers 
contain less than one excitation per channel on average. 
The characteristic measurement time r m will then be on 
the order of the time-of- flight r m ps i/vp of an excitation 
pair through the sample, where t is the average path 
length of the interferometers and vp is the Fermi veloc- 
ity of the ballistic excitations. An observation time of T 
at the drains Di and D 2 therefore roughly corresponds 
to n fa T jr m individual measurement events. 

The contextual values can be used to provide an up- 
per bound for the number of measurement events for 
a desired root-mean-square (RMS) error in the estima- 
tion of the average. Specifically, to estimate the aver- 
age (03) from a sequence of n random drain absorptions 
{di,d,2,--- ,d n ), where di € {Di,D 2 }, one can use an 
unbiased estimator for the average, 

1 - 

E[**] = ~y>*, (37) 

i 

that is defined in terms of the contextual values assigned 
to each measurement realization. As n — ► 00, the esti- 
mator ([57)) converges to (<r|) = ar) 1 Pp )l + a D2 P D2 . The 
mean squared error (MSE) of this estimator is given by 
the variance of the contextual values over the number of 
measurements, 

a 2 n P D , + a 2 Pjy, - (ctI) 2 

MSE Ecx J = Dl 1 ° 2 2 (38) 

n 

Hence, the RMS error ^/MSE[E] scales as l/^/n and im- 
proves with an increasing number of measurements. 

Without prior knowledge of the state, a reasonable up- 
per bound one can make for the MSE is the norm-squared 
of the contextual values over the number of measure- 
ments, 

ct n + a 2 n 

MSE[E[ct,]] < Dl ° 2 . (39) 

It then follows that to guarantee a maximum desired 
RMS error e one needs an observation time on the or- 



der of, 

a 2 n + a 2 n 
r~r m Dl 3 D \ 40 

As the measurement becomes more ambiguous the con- 
textual values become more amplified and so lengthen 
the observation time necessary to achieve the RMS error 
of e. For a strong measurement the upper bound on the 
observation time is T « 2r m /e 2 . 

IV. CONDITIONED MEASUREMENTS 

To gain further insight into the which-path informa- 
tion, we can condition the measurement on the subse- 
quent absorption of the system excitation at a specific 
system drain. To do this we must obtain the joint trans- 
mission probabilities for pairs of detector and system 
drains. Conditional probabilities can then be defined in 
terms of the joint and single transmission probabilities. 

As pointed out by Kang2& these probabilities are ex- 
perimentally accessible in the low-bias regime through 
the zero-frequency cross-correlation noise power between 
a detector drain D £ {Di,D2} and a system drain 
Se{S 1 ,S 2 }, 

Sd,s « 2^(PsA e f) ~ P S (E F )P D {E F )). (41) 

Hence, knowledge of both the average currents (HJ and 
the noise power (|4Tj) allows the determination of both the 
joint and single transmission probabilities. 

A. Joint Scattering 



We can determine the joint probabilities directly in 
the scattering model by rewriting (|10p in the basis of the 
system drains using ([3]), 




up to the same global phase as in (TTU)) . 

The relevant joint scattering amplitudes are, 



£>i,Si — e - r 1 r 2 r 1 r 2 -|- i 2 r 1 r 2 e 

C DlA = e i( *=+^ \rtrir\tl +t(4r{t 2 : e t ^ d+ ~ (} +r(r 2 l tlr°e^ 3 +tUtt s 1 r s 2 e li 'l ,d+ ' l > 3) 
C D2<Sl = e l ^2 +xi ) \ r f t d 2 rj r s 2 + tf r d 2 r\ r s 2 e l ^"+^ + r d t d 2 t\ t s 2 e^' 



r d r d t\t\e^ 



t d t d t{t s 2 e* 



1> d +4> s ) 



+ 4*2*1%** 



+0°) 



(44a) 
(44b) 
(44c) 
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c 



D 2 ,S 2 



Mi +SI ) 



rf 4 r{ t s 2 + tf rf r{ t s 2 e 1 ^"^ + rf 4 t{ r s 2 e**° + if r d 2 t\ r s 2 e l W d +^ 



(44d) 







FIG. 6. (color online) The joint probability Pdx.Sx (|46a[) with 
system tuning <f> s = as a function of the detector tuning (f> d 
and the which-path information 5f , shown for efficient detec- 
tion V d = 1, balanced system drains e| = 1, and the coupling 
phases 7 = {0, 7r/4, tt/2, tt}. 



FIG. 7. (color online) The joint probability Pd 1 ,s 1 (|46a|) with 
system tuning (f> s = tt/2 as a function of the detector tun- 
ing 4> d and the which-path information 8{ , shown for efficient 
detection V d = 1, balanced system drains e| = 1, and the 
coupling phases 7 = {0, 7r/4, it/2, n}. 



The joint probabilities for absorption in detector drain 
D E {Di,D 2 } and system drain S € {Si,<S2} can then 
be understood as either expectations of joint projections 
\S)(S\ (g) \D)(D\ under the joint state ([43]), or, equiva- 
lently, as expectations of system projections |S')(S'| under 
the measured reduced system state Mu\ip s ), 



P D ,s = |(5,D|4-'")| 2 = \(S\M D \r)\ 2 = \C D ,s\ 2 - (45) 



These joint probabilities can be written explicitly in 
terms of the parameters ((4}, (j23|) . and ([28]) as, 



D u Si = 



p n lt s 2 = 



\ (P+P+ + V d V s A ds (46a) 

_ yd ^ips + + T d (^ + gs )) 

-V s (A s l3 d -T s (Sf + Sf))\ 

\(piP- - V d V s A ds (46b) 
- V d (A d (3i_ + T d (SI ~ 5 S 2 )) 
+ V S (A s p d + -T s (Sf + <*£))), 



Pd 2 , Si = \ (ptF + - V d V s A ds (46c) 
+ V d (A d f3 s + + T d {Sf + S s 2 )) 
-V s (A s /3l-r s (sf -sf))\ 

Pd 2 ,s 2 = J (p*PL + V d V s A ds (46d) 
+ V d (A d pt + T d (Sf - S s 2 )) 

+ ys ( A s f3 d_ T s (gd 

where we have defined the additional parameter A ds as, 



A 



ds 



COS COS < 



rds • 1 „• ; / 

= sin — sm 



ds ^ 



(47a) 
(47b) 
(47c) 



For illustration purposes, we have plotted the joint prob- 
ability Pd 1 ,s 1 (|46a|) for several parameter choices in Fig- 
ure [6] and Figure [7] 

The parameter A ds E [—1,1] represents the joint in- 
terference between the system and detector. The pa- 
rameter T ds E [—1,1] is the portion of the joint inter- 
ference that depends explicitly on the coupling phase 7 
and the difference between the tuning phases 4> ds . As 




FIG. 8. (color online) The conditional probability P Dl ^g 1 — 
Pd 1 ,s 1 /Ps 1 with system tunings (p s — (left) and 4> a = 7r/2 
(right) as a function of the detector tuning (j> d and the which- 
path information 6{, shown for efficient detection V d = 1, 
balanced system drains e§ = 1, and the coupling phases 7 = 

{0,7T}. 

the coupling 7 — » 0, then T ds — > and the joint in- 
terference reduces to a decoupled interference product 
A ds — > cos 4> d cos 4> s . As the coupling 7 —> ir, then 
pds _^ cos 0ds an( j the j i n t interference will be maxi- 
mally coupled A ds — ► — sin <fi d sin <f> s . 

We can marginalize the joint probabilities (pQj)) to 
obtain both the detector probabilities (|27|) as Pd = 
Ss-^d.s an d the system probabilities (|2T))) as Ps = 
Pd.s- Furthermore, we can construct the conditional 
probabilities Pd\s f° r absorption in a detector drain D 
given an absorption in a system drain S, as well as the 
conditional probabilities Ps\d f° r absorption in a system 
drain S given an absorption in a detector drain D, 

Pd\s = (48a) 
Ps 

Psp = (48b) 
Pd 

For comparison with the joint probabilities, we have plot- 
ted the conditional detector probability Pd^s^ in Fig- 
ure [S] for several parameter choices. 



B. Quantum Erasure 

We can use the conditional system probabilities Ps\d 
to clarify the phenomenon of quantum erasure^'— 
which has also been explored in this system by Kang 38 . 
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FIG. 9. (color online) Quantum erasure. Any system in- 
terference in the system drain probability Ps 1 (upper-left, 
black) is completely destroyed at strong coupling 7 = n (gray, 
dashed) . Conditioning on the detector drains recovers phase- 
shifted interference in the conditional probabilities Ps 1 \n 1 and 
Psi|d 2 (bottom), but with a visibility that is dependent on 
the measurement ambiguity as controlled by the detector tun- 
ing phase <j> . The contextual values {ao 1 ,aD 2 } (upper-right) 
diverge for maximum ambiguity and reduce to the eigenval- 
ues of the which-path operator 03 for zero ambiguity. The 
plots are shown for efficient detection V = 1, strong cou- 
pling 7 = 7T, and maximum system visibility V s = 1. 



Any wave-like interference patterns in the system drain 
probabilities will degrade with the coupling phase 7, as 
shown in the upper-left of Figure [9] for Pg t and maxi- 
mum system visibility V s — 1. At strong coupling 7 = tt 
the wave-like interference will be completely destroyed. 
However, parts of the interference may be restored by 
conditioning the drain on an appropriate detector mea- 
surement. 

To restore the interference in the system statistics, the 
detector must make an ambiguous measurement, as we 
shall see. Strong coupling destroys the interference in the 
unconditioned system statistics by recording the which- 
path information in the detector state via the coupling 
phase 7. As the which-path information is available in 
the detector state for later collection, at least in principle, 
the total reduced system statistics must reflect the degree 
of potential information acquisition. However, such infor- 
mation in the detector state has not yet been extracted 
classically since the detector drains have not yet been 
probed; hence, the information in the state only indicates 
the potential for later extraction of which-path informa- 
tion at the detector drains. A partially ambiguous mea- 
surement extracts some of that potential information and 
erases the rest; a completely ambiguous measurement ex- 
tracts no information and erases all of the potential for 
later extraction in the process. The recovered interfer- 
ence in the conditioned system statistics reflects the era- 
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sure of the information acquisition potential by the am- 
biguous measurement, even though the total statistics of 
the reduced system are unchanged by probing the detec- 
tor. 

As discussed in Section ITlI El the detector phase <fi d de- 
termines the ambiguity of the measurement under such 
strong coupling, so the degree of possible erasure will 
also depend on the detector phase. We can see the de- 
pendence of the interference recovery on tfi d in Figure [9] 
in the lower two plots. As the detector phase cf> d varies 
from to 2ir, the conditional probabilities -PsilDi an d 
Psi\d 2 continuously vary from flat particle-like statistics 
to complementary wave-like interference patterns. The 
visibilities of the complementary interference patterns di- 
rectly depend on the measurement ambiguity, as can be 
seen in the plot of the contextual values in the upper- 
right of Figure [9] Maximum visibility corresponds to 
maximum measurement ambiguity where the contextual 
values diverge; zero visibility corresponds to zero mea- 
surement ambiguity where the contextual values reduce 
to the eigenvalues of the which-path operator. We also 
note that the effect of the additional coupling evolution 
pOep creates a ir/2 phase shift in the interference pattern 
that cannot be erased since it is not part of the informa- 
tion extraction of the measurement. 

Such erasure under strong coupling has an intuitive 
analogy to an optical double-slit experiment (that you 
can even try at home!) 4 ^ as shown in Figure [TO] In 
the optical equivalent, a coherent beam of light passes 
through two slits and displays an interference pattern on 
a remote screen. However, if an experimenter tags each 
slit with horizontal and vertical polarizing filters, then 
the which-path information of the light can later be ex- 
tracted from the polarization degree of freedom; hence, 
the total interference pattern on the remote screen will be 
destroyed. The experimenter can subsequently condition 
the measurement on a particular polarization by plac- 
ing another polarizer after the two slits oriented at some 
angle 9 relative to horizontal. If the conditioning polar- 
izer is oriented horizontally or vertically, then the path 
measurement will be unambiguous and extract all which- 
path information, so no conditioned system interference 
will be recovered. However, if the conditioning polarizer 
is oriented diagonally, then the path measurement will 
be completely ambiguous and the potential which-path 
information will be erased, recovering all of the inter- 
ference in the conditioned statistics. In the electronic 
version, the system excitation plays the role of the light 
beam, the relative coupling phase 7 records the potential 
path information, the conditioning on a particular detec- 
tor drain plays the role of the polarizer, and the tuning 
phase (f> d selects the conditioning basis analogously to the 
angle controlling the ambiguity of the measurement 
and degree of erasure. 

Furthermore, one could in principle implement a 
delayed-choic o 42 i 43 ' 45 version of the quantum erasure by 
placing the detector drains much further away in the sam- 
ple than the system drains. The interaction phase 7 could 



a) 



b) 




FIG. 10. (color online) Analogy to optical quantum erasure of 
two-slit interference, a) A coherent light beam passes through 
two slits tagged with horizontal (h) and vertical (v) polariza- 
tion. After conditioning on a subsequent polarizer oriented 
at an angle 9 to horizontal, an interference pattern can be re- 
covered with a visibility that depends upon 9. b) A coherent 
electron passes through a QPC and is tagged by a detector 
electron via an interaction phase shift 7. After conditioning 
on the drain D2 after a subsequent detector QPC that forms 
an MZI with tuning phase <f> d , an interference pattern can be 
recovered in the system drain Si as the system tuning 4>" is 
varied with a visibility that depends upon <f> d . 



be recorded and the system excitations collected, upon 
which a controlled change in the magnetic field could set 
the tuning phase of the detector. Upon conditioning the 
data, the interference would reappear according to which 
tuning phase had been chosen after the system excitation 
had already been collected. 

We stress that the erasure of the potential which-path 
information and recovery of the system interference will 
be apparent only when conditioning the collected data. 
Without conditioning, even completely ambiguous mea- 
surements under strong coupling will destroy the system 
interference. The interference patterns recovered by con- 
ditioning on the detector drains will be complementary 
to each other in such a case and thus cancel in the un- 
conditioned statistics. 



C. Conditioned Averages 

We can also use the conditional probabilities to con- 
dition the averages of the which-path measurement on 
a subsequent system drain absorption. To do this, we 
weight the conditional detector probabilities (|48p with 
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the contextual values for the measurement (l25l) . 55 

s (al) =J2 a D P D\s- ( 49 ) 

D 

In terms of the various characterization parameters, these 
take the explicit form, 

sA°l) - {si +51 V s ^Pj , (50a) 

52 ( ^ } = 2pJ 2 { 5{ ~ 51 + yS T^ ) ' (50b) 

where there is a joint interference contribution to the 
conditioned average, 

E ds = A ds -(A s -SfT s ) (51) 

The joint interference simplifies considerably in the 
special case of efficient detection V = 1, 

=J^2sin^cot(!+0 d )cos(^-^). (52) 

This case is plotted in Figure [IT] 

The interference term is scaled by the visibility of the 
system interference V s , which measures the wave-like be- 
havior of the system excitation. Any wave-like behavior 
of the system leads to a contribution to the conditioned 
averages that depends on properties of the correlated de- 
tector as well, due to the joint interference. Hence, the 
which-path information of the system cannot be decoupled 
from the detector that is measuring it in general. 

The conditioned averages properly obey the consis- 
tency relation (<r|) = ^2 S s(<r|) Ps, and are bounded by 
the contextual values (|25[) . Since the contextual values 
are usually larger than the eigenvalues of <r| due to am- 
plification from measurement ambiguity, the conditioned 
averages can counter-intuitively lie outside the eigenvalue 
range. In the weak coupling limit 7 — > 0, such con- 
ditioned averages can become detector-independent and 
converge to weak valuesf 55 ^ 63 as we will show later. 

However, for any macroscopic property such condi- 
tioned averages will always lie inside the eigenvalue range, 
even when measured ambiguously. Specifically, the eigen- 
value range constraint for conditioned averages has been 
shown to be equivalent to a generalized Leggett-Garg 
constrain t 56 ! 72 ' 73 that must be satisfied for any non- 
invasively measured, realistic property. As such, any vi- 
olation of the eigenvalue range in a conditioned average 
can be seen as a signature of nonclassical behavior stem- 
ming from quantum interference. 

D. Deterministic Measurement 

If either of the system QPCs is fully transmissive or 
reflective then V s = 0, the system interference vanishes, 




FIG. 11. (color online) The conditioned average s i (o'I) (|50a[) 
as a function of the coupling phase 7 and the which-path 
information 5( , shown for efficient detection V d = 1, balanced 
system drains e| = 1, and tuning phases (j> d , <j) s — {0, 7r/2}. 

and the conditioned averages (|50l) reduce to ±1 for any 
coupling strength. In such a case the excitation path 
is deterministic and the system displays purely particle- 
like behavior. The post-selection perfectly determines 
the path, and the which-path measurement made by the 
detector will always agree with the post-selected value. 

E. Strong Coupling 

For the case of strong coupling 7 = 7r and an efficient 
detector V d = 1, the conditioned averages ([ST))) reduce 
to, 

s, (^3> "> ^-J^ + ^ tan <p d sin <f , (53a) 
P+ P+ 

s^M^^-^tan^sin^. (53b) 

The tuning phase <fi d is the sole detector parameter that 
specifies the ambiguity of the measurement. 

If the measurement is also unambiguous <fi d — > nir, 
then the interference contribution vanishes. The condi- 
tioned averages become the detector-independent quanti- 
ties Sl (a s 3 ) -> (St +S s 2 )/f3 s + and s 2 <a|> -> (Sf -5f)//3I 
that always lie in the eigenvalue range. A strong which- 
path measurement made by the detector therefore forces 
the system excitation to display particle-like conditioned 
statistics. 

However, even with strong coupling any ambiguity in- 
troduced in the detector will lead to quantum erasure 
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FIG. 12. (color online) The weak limit 7 — > of the con- 
ditioned averages Si (03) (top) and s 2 (o"I) (bottom) for the 
distinct cases of weak measurement <f> d 7^ (left) given in (|54[) 
and semi-weak measurement d = (right) given in (I56|l as a 
function of the system tuning cf> a and the which-path informa- 
tion 8{. The values are shown for efficient detection V = 1 
and balanced system drains e| = 1. 



that recovers interference in the conditioned statistics 
of the system, as discussed in Section IIVBI Such recov- 
ered interference can take the conditioned averages of the 
which-path information outside the eigenvalue range. For 
an almost completely ambiguous measurement the tun- 
ing phase will deviate from ir/2 only by a small angle 
5<f> d . Since tan(7r/2 + 5<j> d ) = -1/ '5<f> d + 0{8(f) d ), the inter- 
ference contribution will dominate, and the conditioned 
averages will diverge. 



F. Weak Coupling Limit 

The weak coupling limit 7 — » of an efficient detector 
V d = 1 leads to conditioned averages (j5Tj)) that generally 
depend on the detector tuning <fi d , as anticipated during 
the discussion in Section flll Fl For the weak measurement 
case when the detector tuning is not an integer multiple 
of 7T, then as 7 — > the joint interference term in the 
numerators of ([50)) vanishes, sin(7/2) cot(7/2 + <p d ) — > 0, 
yielding the detector-independent quantities, 



P s , - V s cos ( 



SI 



P s _ + V s cos< 



(54a) 
(54b) 



In contrast with the unambiguous case, system interfer- 
ence remains in the denominators. 



These expressions match the real parts of the weak 
value expressions defined in Aharonov et ali^ 



S 1 \0- 3 )w - 



V s sin <j) 8 



(3 S + - V s cos < 



(55a) 



(55b) 



V s sin < 



(Sim ' 

Si +5 S 2 
/3 S + - V s cos <fr 

si - sj 

P s + + V s cos <j) s ^ 1 P s + + V s cos <f) s 



As pointed out in Dressel et al<^ any additional unitary 
coupling evolution such as (|30ej) could in principle af- 
fect the convergence of the conditioned averages to these 
detector-independent weak value expressions. However, 
in this case the limit is unaffected and the standard ex- 
pressions are recovered. 

For the semi-weak measurement case when <\> d = 
nn with n an integer, then as 7 — > the joint in- 
terference term in the numerators of (|50p reduces to 
sin(7/2) cot(7/2 + nn) — > (—1)"; therefore, system in- 
terference remains in both the numerators and denomi- 
nators, yielding the modified expressions, 



Si (03) 



Sf +51 -(~l)"F s cos(f 

^-V'cos^ ' [ j 

Sf - j| +(-l)T'cos^' 

p s _ +V s coscj) s 1 ' 



The integer n selects which detector drain is projective. 
Hence, unlike the weak values (f5"5"j). these "semi- weak 
values" do not conform to general detector-independent 
expressions, but explicitly depend on the details of the 
measurement operators (|30|) . 

We emphasize that the detector tunings 4> d = nir are 
critical points around which the 7—^0 limit is unstable, 
meaning that any laboratory approach to the weak cou- 
pling limit of the conditioned averages ([517)) can approxi- 
mate either the weak values (|5~4"|) or the semi-weak values 
([55)) depending on the relative magnitudes of 7 and (j) d . 
Hence, the limiting values will compete with each other 
as 7 becomes small and (j) d approaches a critical point. 
The difference between these limiting cases is plotted in 
Figure [T2l using the choice n — 0. 



V. CONCLUSION 

This work describes in detail how one electronic Mach- 
Zchndcr interferometer may be used as a detector for 
the which-path information of a second Mach-Zchnder 
interferometer. We use a combination of the Landauer- 
Buttiker formalism of coherent charge transport and the 
contextual values formalism of generalized measurement 
to demonstrate that inducing an interaction phase shift 
in the joint state for adjacent arms of the interferometers 
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allows only the which-path information to be inferred. 
For clarity of discussion, we use a simplified two-particle 
model of a Coulomb-induced joint phase shift in the low- 
bias regime to couple the interferometers and perform the 
measurement; however, the technique may be extended 
to a microscopic many-body model of the coupling inter- 
action as well. The efficiency of the which-path mea- 
surement depends on the interference visibility in the 
detector, demonstrating wave-particle complementarity 
between the detector and the system being measured. 
Moreover, we find that the total ambiguity of the mea- 
surement depends not only on the coupling strength and 
the measurement efficiency, but also on the tuning phase 
of the detector interferometer. Such additional measure- 
ment ambiguity leads to quantum erasure of the poten- 
tially extractable which-path information in the detector 
state; the erasure can be exposed by conditioning the 
system statistics on specific detector outcomes. 

The contextual values, or generalized eigenvalues, of 
the which-path operator become amplified from the 
eigenvalues to compensate for any ambiguity in the mea- 
surement. The analytic form of the amplification gives 
direct insight into the sources of the ambiguity and in- 
dicates the proper background removal for the detector. 
The ambiguity compensation allows correct which-path 
information to be obtained on average, despite ineffi- 
ciency of the measurement. The contextual values also 
allow the principled construction of conditioned which- 
path averages for any coupling strength. In general, such 
conditioned averages contain information about both the 
detector and the system being measured due to the joint 
quantum interference, but in the strong or weak coupling 
limits they can converge to detector-independent quanti- 
ties. In the weak coupling limit, we can obtain not only 
detector-independent weak values from the conditioned 
averages, but also detector-dependent semi-weak values 
at critical points of the detector tuning. The two distinct 
limits will compete depending on the relative magnitudes 
of the coupling strength and the detector tuning phase. 

In addition to providing a simple but nontrivial theo- 
retical example of a quantum mechanical detector inter- 
acting with an identical quantum mechanical system, the 
coupled electronic interferometers that we consider can 
be realized with current technology. Furthermore, our 
results will hold with only minor modifications for an op- 
tical equivalent using a two-photon interaction (such as 
a nonlinear crystal) to induce the relative phase shift be- 
tween copropagating interferometer arms. Since the in- 
formation being measured can be easily controlled by the 
experimenter in both situations, our measurement proce- 
dure can be verified experimentally. The same techniques 
used here can be used to calibrate and analyze quantum 
mechanical detectors for less controllable system infor- 
mation in future experiments. 

We acknowledge helpful discussions with I. Ncder and 
support from the National Science Foundation under 
Grant No. DMR-0844899, and the US Army Research 
Office under grant Grant No. W911NF-09-0-01417. 



Appendix A: Interaction Phase 

It may not be apparent that an extended Coulomb in- 
teraction between scattering excitations in adjacent edge 
states can result in an additional joint phase accumula- 
tion 7 without destroying the phase coherence. Indeed, 
in momentum-space the Coulomb interaction explicitly 
involves energy exchange between the adjacent excita- 
tions, which would seem to imply that phase disruption 
would occur. To assuage such concerns, we shall solve a 
simple model of Coulomb edge state interaction for the 
two-excitation amplitude in the co-propagating region. 
We shall see that it is sufficient to keep the total energy 
of both excitations constant in order to obtain a joint 
phase accumulation over the interaction length; the fact 
that the excitations may exchange energy between them 
does not disrupt the joint phase coherence. 

Consider the longitudinal part of a two-particle am- 
plitude ip(xi,X2,t) that describes chiral copropagation 
at a speed v along linear edge channels separated by a 
fixed distance d. In the absence of any Coulomb inter- 
action between the channels, tp must satisfy the effective 
Schrodingcr equation, 



ihd t ip = -^(<9xi + d X2 )V>. 



(Al) 



For a fixed total energy E, the general solution of this 
equation has the form, 

ip E {xi,x 2 ,t) = i>o{t)tpi(xi,x 2 )x{x2 - xi), (A2a) 

Mt) = iPoexp(-ij^j , (A2b) 

ipi(xi,x 2 ) = exp (jj^&i +x 2 ) S j , (A2c) 

where \ is an arbitrary function of the difference of the 
coordinates. 

By choosing an initial boundary condition to be a prod- 
uct state of excitation scattering states at distinct ener- 
gies, 



E±=E(l±p), 



(A3) 



where p G [—1,1] such that the energy matching con- 
dition E+ + E~ = 2E is satisfied and E is correctly 
quantized, we fix \ to find the general product form for 
a fixed joint energy of IE, 

il>p(2E,x 1 ,X2,t)=t(E+,x 1 ,t)Z(E-,x 2 ,t), (A4a) 

/ E 2r \ 
Z(E,x,t)=Z e X pl-ij;{t-—)). (A4b) 

As expected, the channel amplitudes are completely de- 
coupled in the absence of interaction and independently 
phase-coherent. 

For a low-biased source, each single-particle energy will 
be approximately the Fermi energy _E+ w E~ w Ep, and 
the propagation speed will be the Fermi velocity v = vp. 
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Therefore, in the absence of interaction each particle will 
accumulate a dynamical phase, 



cj>(E F ,L) = ^L, 



(A5) 



after a propagation length L, leading to a total joint dy- 
namical phase of AEpL/hvp- 

If the co-propagating excitations are instead allowed to 
interact via a screened Coulomb potential, the effective 
Schrodinger equation must be modified to, 



ihdtip 



ihv 
~2 



(d Xl +d X2 )i> + —e- r ' x iP, (A6a) 



r = \J d 2 4- \x 2 - £i| 2 , 



(A6b) 



where A is the screening length, r is the interaction dis- 
tance that depends in the interacting region on the differ- 
ence \x 2 — X\ | between the coordinates as well as the dis- 
tance d between the edge channels, and a is the Coulomb 
interaction constant in appropriate units. This linear 
equation decouples in the coordinates yi = x± + x 2 and 
j/2 = X 2 — x\, so it may still be exactly solved. 

For a fixed joint energy E, the general solution has the 
form, 

i>E{x\,x 2 ,t) = ipo{t)^i{x u x 2 )x{x2 - xi), (A7a) 
i>i{xi,x 2 ) = exp (ik(E, x±, x 2 )(xi + x 2 )j , (A7b) 

(A7c) 



1 



k(E,x 1 ,x 2 ) = -^(E-— e~ r ' x ) 
nv r 



where x/jq is the same as in (|A2I) and x 1S another arbitrary 
function of the difference of coordinates. The Coulomb 
potential thus gives an effective position-dependent shift 
to the joint wave-number k for the amplitude, which will 
affect the dynamical phase accumulation for the joint am- 
plitude. 

If we demand that for r A the general solution (IA7|) 
should reduce to the noninteracting solution (|A4[) . then 
we find the simplest form, 



i> p (x 1 ,x 2 ,t)=iP P (2E,x 1 ,x 2 ,t)e- i ^' x *\ (A8a) 

(A8b) 



ae" _ r/x xi + x 2 

7(xi,a;2) = — e 1 . 

fir v 



where tfj p is the noninteracting solution (|A4|) . The net 
effect of the Coulomb interaction between the channels 
is thus to contribute a position-dependent (but energy- 
and time- independent) phase 7(2:1, x 2 ) that entangles the 
coordinates of the channels. Any remaining correction 
factor x{x 2 —X\) to this simple form must satisfy x(0) = 
1, so we will safely neglect it in what follows. 

The excitations are collected at ohmic drains at fixed 
positions x\ = L\ and x 2 = L 2 of the coordinates, so 
the detected phase j(Li,L 2 ) will be fixed by the geom- 
etry and stable for any pair detection. Scattering states 
with fixed energy such as (|A8j) are stationary and ex- 
tended throughout the interaction region, which explains 



the geometric nature of the interaction-induced phase. 
The square of the wave-function \ip p (2E, L\, L 2 , t)\ 2 in- 
dicates the (typically small) probability that the excita- 
tions will be detected simultaneously at any particular 
t at the drains. However, by using a coincidence post- 
selection or by engineering a correlated initial scattering 
state, one can in principle restrict the bulk of the mea- 
sured detections to be coincident. 

If we further assume that before and after a co- 
propagation length L where the channels are a fixed dis- 
tance d apart both edge channels rapidly split away from 
one another, then to good approximation the Coulomb 
interaction only affects the region of length L. After the 
interaction region, the equation of motion for the state 
will effectively revert to (|Al|) . restoring noninteracting 
dynamical phase accumulation similar to (|A5[) . Hence, 
the amplitude for jointly detecting the excitations at 
ohmic drain positions L\ > L and L 2 > L will contain an 
extra dynamical phase that is accumulated by the joint 
state only within the interaction region, 



ipp {Li ,L 2 ,t) =ip p (2E,Li,L 2 ,t) e _n , 
, r t \ «e 2 ai\2L 
nd v 



(A9a) 
(A9b) 



Moreover, the joint state may be further scattered after 
accumulating the joint interaction phase but before the 
joint detection, as indicated in the main text. We see that 
for such simultaneous detection the interaction phase 7 
will be linear in the interaction length L and therefore 
should be tunable in principle. 

If the detections are not simultaneous, then one exci- 
tation will be detected in a drain at a time t\ , followed by 
the second excitation at a later time t 2 . The joint state 
will therefore be successively collapsed by the drain de- 
tections, which will introduce an additional relative phase 
factor due to the discrepancy in detection time. Specifi- 
cally, for xi, x 2 > L the joint state will have the form, 



4> p (x 1 ,x 2 ,t) = i/j p {2E,xi,x 2l t)< 



(A10) 



with the accumulated interaction phase 7 as in (|A9|) . De- 
tection of the first excitation at L\ at t\ then collapses 
the joint state to, 



ip p (2E,Li,x 2 ,t) 



-*7 



P{L u ti) = I dx 2 \$ p {L 1 ,x 27 t 1 )\ 2 . 



(Alia) 
(Allb) 



Evolving the remaining single particle state to time t 2 
and then detecting the second excitation at L 2 produces 
the amplitude, 



7 „ ^ P (2S,L 1 ,L 2 ,t 1 ) 
Wv = : — e 



-i(7+B(t 2 -ti)/fi) 



(A12) 



The only difference between the sequential detection am- 
plitudes and the joint detection amplitude (|A9I) is the 
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extra temporal phase that is accumulated between the 
two detections. Notably, the extra temporal phase factor 
appears as a global phase that should not affect the final 
statistics, in contrast to the geometric interaction phase 
7, which can be exposed by further scattering before the 



sequential detection as in the main text. Allowing for 
fluctuations in 7 as in ^III Gl will account for geometric 
uncertainty in the interaction length, as well as the pair 
injection frequency of the source. 
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